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Synopsis 

Lecture 1 : 

• Review of trigonometric identities 

• Fourier Series 

• Analysing the square wave 

Lecture 2: 

• The Fourier Transform 

• Transforms of some common functions 

Lecture 3: 

Applications in chemistry 

• FTIR 

• Crystallography 
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Introduction 


Chemistry often involves the measurement of properties which are the aggregate 
of many fundamental processes. A variety of techniques have been developed for 
extracting information about these underlying processes. Fourier analysis is one 
of the most important and is very widely used - eg: in crystallography. X-ray 
adsorbtion spectroscopy, NMR, vibrational spectroscopy (FTIR) etc.. As it involves 
decomposition of functions into partial waves it also appears in many quantum 
mechanical calculations. 


A Little Trigonometry 

You will need to be able to manipulate sin() and cos() in order to understand 
Fourier analysis - a good understanding of the UK’s A-level Pure Maths syllabus is 
sufficient. Here is a brief reminder of some important properties. 

Units angles are typically measured in radians: 0 — 360°is equivalent to 0 — 27 t 
radians 

Cos and sin curves look like this: 
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Both sin(x) and cos(x) are periodic on the interval 2tt and integrate to 0 over a full 
period, ie: 

/* + 7T /*+7T 


/ cos (x) dx = 


sin (x) dx = 0 


Wavelength 

It, should be clear that sin(2x) repeats on the interval 0 —> tt and sin(3x) on the 
interval 0 —> 27r/3 etc. In general sin(nx) and cos(nx) repeat on the interval 0 —> 
27r/n. The repeat distance is the wavelength A and so in general, \ = 2Tr/n. 

The discrete family of functions sin(nx), cos(nx) are all said to be commensurate 
with the period 27r- that is, they all have wavelengths which divide exactly into 27r. 

The function sin(kx) for some real number k has an arbitrary wavelength A = 27 t /k. 

k is usually referred to as the wavevector. 

Note: A simple Mathermatica notebook, trig_l.nb, is provided with the course 
and can be used to play with sin and cos functions. 


Fourier Series 

The idea of a Fourier series is that any (reasonable) function, f(x), that is peri¬ 
odic on the interval 27 t (ie: f(x + 27m) = f(x) for all n) can be decomposed into 
contributions from sin(nx) and cos(nx). 

The general Fourier series may be written as: 


f(x) = -^ + a\ cos (x) + a 2 cos ( 2x ) + <23 cos (3a;) + ... + a n cos (nx) 

+ b\ sin (a;) + 62 sin (2a;) + 63 sin (3a;) + ... + b n sin (nx) (1) 

Note: 

1 . cos (nx) and sin (nx) are periodic on the interval 27 t for any integer n. 

2. The a„and b n coefficients measure the strength of contribution from each 
“harmonic”. 

Orthogonality 

The functions cos (nx) and sin (nx) can be used in this way because they satisfy the 
following orthogonality conditions: 


3 



f + TT 
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cos (mi) sin (nx) dx 
cos ( mx ) cos (nx) dx 


(•+7T 


sin (mx) sin (nx) dx 


0 for all to, n 

0 m=f n 

2i r to = n = 0 

7r (to- = n) > 0 

0 to ^ n 

7r (to- = n) > 0 


Note that the integrals only need to extend from — 7r to +7r (or any other period of 
27t) as the functions simply repeat outside this range. 

These conditions can be proved quite readily but it is relatively easy to see why 
they are true graphically. 


• cos(mx) sin(nx) ? This is obvious (!) if you plot cos(x) and sin(x). 

• sin(mx) sin(nx) n to? It is easier to see why it is true by picking a special 
case; say the integral of sin(2x) sin(x) and plotting. 


Sin^x} x 3in(£x) 
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The symmetry of the plot makes it clear that an integral of this function over any 
period of 27r will yield 0. 


• Is it obvious that this will be true for all cases when n ^ ml 

• Also for the case cos(mx)cos(nx) ? 

• Can you prove it in the general case ? 


Note: A simple Mathermatica notebook, trig_l.nb, is provided with the course 
and can be used to play with these products and integrals. 


cos(mx) cos(mx) - ie: the case when n = 


/ + 7T /*+7T 

cos 2 (mx) dx= - (1 + 2 cos (mx)) 

-IT J — 7T " 


]dx= 2 


sin (2 mx) 


2 TO 
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Finding the coefficients 


As was shown in the lecture the orthogonality conditions allow us to pick off values 
for all of the coefficients. Multiplying the whole Fourier series by l,cos(na;)or 
sin (nx) and integrating over a complete period leads to terms which are zero apart 
from one which corresponds to the coefficient ao,a n or b n respectively, that is: 


1 /'+ x 

a 0 = - f(x)dx 

1 f +n 

a n = — / f (x) cos (nx) dx 

K J- 7T 

1 f + 7T 

b n = — / f (x) sin (nx) dx 

If f(x) is “well behaved” we can perform these integrals and obtain the Fourier 
decomposition of f(x). 

Note: Well behaved in this context means that the function obeys the Dirichlet 
conditions. 


1 0 < X < 7T 

0 —7T < x < 0 

f(x + 2tt) 
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An Example 

Consider the square wave: 

f(x) = 
f(x) = 


This appears to be a difficult case - the rather angular square wave does not look 
as if it will be readily expanded in terms of sine and cosine functions. 

The coefficients in the expansion can be determined from the formulae given above. 

cio is determined by: 
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2 /*+7T ^ /*+7T 

ao = — / f(x)dx = — 1 dx = — X7r = l 

TT ^— 7T 7T 7o ^ 

where the restriction of the integral to the region 0 —> tt is simply because /(x) is 
zero in the region —tt —> 0. 

Similarly: 


1 /' +7r 

a n = — / lx cos (nx) dx = 0 
^ Jo 


just draw cos (x)to see that its integral from 0 —» tt is zero, so 


a n = 0 for all n 


For the b coefficients we have, 


C+n 


bn — 


1 x sin (nx) dx 


— cos (nx) 
n 


J o 


(1 — cos (wr)) 


but, 


cos (rwr) = +1 n — even 

= —1 n — odd 


so, 


b n = 0 n — even 


having determined all of the coefficients we can write the series for /(x) as: 

... 1 2 / . sin (3x) sin (5x) \ 

/(*) = - + -(*,„ W + -i-U-LJ + ...) 

The sum continues to an infinite number of terms. We can see how it converges to 
the square wave by plotting the truncated sum containing a finite number of terms 
- lets call the sum containing n-trigonometric terms /„(x) then /o (x) = jps just the 
average value of the square wave. 

/i(x) = | ^ sin (x) is plotted below 
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which is the “best” approximation that can be made using just a constant and a 
sine wave - not great. 

fz(x) = \ + ^ sin ( x ) + j^sin(3x) looks like this: 



which is beginning to look more like a square well. 

The weight in each contribution is falling and with each additional term the fine 
detail of the square wave is being refined. 

Including 50 terms, / 50 (a;), we get 





which is a pretty decent approximation to the original square wave. 

Note that the little spikes at the edge of the square wave are present even after 
including many hundreds of terms (although they become finer and finer) they are 
a consequence of trying to describe a discontinuous step function with smooth sine 
waves - this was noticed and studied by the mathematician JW Gibbs in the late 
1890’s. 

A more compact notation 

In many applications you will find that a more compact notation is used for the 
Fourier series. Using the identity 


e ld = cos ( 9 ) + i sin ( 9) 

We can write; 

cos (6») = ^ (e l9 + e~ w ) 

sin {0) = ^(e td -e~ ie ) 

Using these relations we can rewrite the Fourier series, equation 1, in the more 
compact exponential notation; 


f(x) 


which can be rearranged as; 


/O) 


y , C n e 
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In the exponential notation the orthogonality conditions are; 


r+TT 


if m = n 


and so the coefficients given by 


r +TT 


2ir 


c f{x)dx 


The original and compact notations are equivalent and the c n coefficients are there¬ 
fore directly related to a n and b n 


1 

Co — 2°° 

Cn = ^ ^ n ) 

C—n = ^ (® n 

This more compact notation is used in almost all applications. 


Fourier Transforms 

Functions of arbitrary periodicity 

The discussion of Fourier Series above dealt with functions periodic on the interval 
27 t (ie: f(x + 2nn) = f(x) for all n). This can be generalised to functions periodic 
on any interval. 

Functions with a periodicity of 2L (ie: f(x + 2 Ln) = /( x) for all n) can be decom¬ 
posed into contributions from from sin (n^)and cos (n^) which are periodic on 
the period 2L. 

The Fourier series may then be written as: 


x a o , firx\ 

f(x)= y + a i cos ) 


+ bi sin 


( 2 


+ o 2 cos (2^) + a 3 cos ^3^ 


L ) 


L 


( t:x\ 
■ ■ ■ + a n cos y n ~j- J 

1 . ( TTX\ 

.. + b n sin [n—J 


( T )+6 2S m(2 T )+i. 3S m(3 T ) 


or 


OO __ 

f(x) = y + ^2 (an cos (npj + 6„sin 

n— 1 

or, in the exponential notation, 

OO 

f( X ) = J2 c ^ inT 

71 — — OO 
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1 

2 L 


+L 


~.f{x)da 


Note: The limits of integration cover a single period of the function which is not 
2L rather than 27 t. 

This allows a function of arbitrary period to be analysed. 


Nonperiodic functions 

Fourier series are applicable only to periodic functions but non-periodic functions 
can also be decomposed into Fourier components - this process is called a Fourier 
Transform. 

Imagine a function that is of a finite extent that is much less than the periodicity, 
2L, as pictured below, 



-L 



L 




If L becomes very large (tends to infinity) then we will have an isolated, aperiodic, 
function. We will use this limiting process to develop the equations for the Fourier 
Transform from the Fourier Series. 

Consider the Fourier Series for this function; 


OO 

f(x) = 

n =—oo 


Consider the limit in which L becomes very large. 
If we define; 


k 


n 


717T 

IT 


then 


OO 

f(x) = Cnelk 

IT— — OO 


and it is clear that for very large L the sum contains a very large number of waves 
with wavevector fc„and that each succesive wave differs from the last by a tiny 
change in wavevector (or if you prefer, wavelength), 


A k = k n+ 1 - k n 


IT 

L 
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As was shown in the lecture, in the limit of large L k becomes a continuous vari¬ 
able, the discrete coefficients, c n , become a continuous function of k, c(fc)and the 
summation can be replaced by an integral and, 


f(x) 

c{k) 


1 r +oo 

— / c(k)e lkx dk 

27T 7-00 


r*+oo 


f(x) 


—ikx 


dx 


These pair of equations are very often rescaled by substituting c(k ) = \Z2nc(k) to 
obtain; 


i r+oo 

f(x) = -= / c(k)e ikx dk 

V Z7T J — oo 

1 /* + °° 

c(k) = -= / f(x)e~ ikx dx 

V Z7T J — oo 

The functions f and c are called a Fourier transform pair - c is the Fourier transform 
of f and f is the (inverse) transform of c. 
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